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We extend previous work showing that violation of the null energy condition implies instability 
in a broad class of models, including gauge theories with scalar and fermionic matter as well as 
any perfect fluid. Simple examples are given to illustrate these results. The role of causality in 
our results is discussed. Finally, we extend the fluid results to more general systems in thermal 
equilibrium. When applied to the dark energy, our results imply that w = p/p is unlikely to be less 
than — f. 



I. INTRODUCTION 

Every spacetime is a possible solution to Ein- 
stein's equations for some particular choice of energy- 
momentum tensor T^ v . It is therefore important to un- 
derstand what general restrictions can be placed on the 
energy-momentum of a physical system. Such restric- 
tions are referred to as energy conditions, and play an 
important role in general relativity. In an earlier paper 
, two of the present authors demonstrated a direct con- 
nection between stability and the null energy condition 
(NEC) 0, T llv n^n v > for any null vector rt l (satis- 
fying Q^ v n il n v = 0). The main results of the study in 
Ref. m are: (1) classical solutions of both minimally and 
non-minimally coupled scalar-gauge models which vio- 
late the NEC are unstable, (2) a quantum state in which 
the expectation of the energy-momentum tensor violates 
the NEC cannot be the ground state (includes models 
with fermions), (3) perfect fluids which violate the NEC 
are unstable. These results suggest that violations of 
the NEC in physically interesting cases are likely to be 
ephemeral. 

Using the Raychaudhuri equation for the expansion of 
a hypersurface null congruence, one can show that null 
geodesies have non-increasing expansion as long as the 
NEC is obeyed Q- This result plays an important role 
in the classical singularity theorems of general relativity 
and in proposed covariant entropy bounds Q| . It also 
implies that NEC violating matter is necessary for the 
construction of Lorentzian wormholes [f| . 

This work is an extension of Ref. . We give some ex- 
plicit examples of the classical stability analysis (Sec. ITT)l . 
additional details regarding the local instability in quan- 
tum states which violate the NEC (Sec. IIH0 . a simplified 
and more general proof of the result that the fermionic ef- 
fective action does not lead to violation of the NEC (Sec. 
II V|) . and a generalization of our analysis of fluids to any 
system in thermal equilibrium (Sec. EJ. In Ref. ||| 
it was shown that models with superluminal modes can 
evade our results in some cases. However, in Ref. it 
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was noted that such models are acausal. In the appendix, 
we show that in causal models, the NEC is sufficient to 
determine stability. In Sec. [H]it is shown that considera- 
tions involving causality or superluminality are necessary 
only in cases where the classical background violates all 
rotational symmetries. Cosmological models of the dark 
energy must exhibit isotropy and hence are covered by 
the earlier results of Ref. . 

The NEC has direct implications for the dark energy 
equation of state, often given in terms of w — p/p. Dark 
energy has positive energy density p and the energy- 
momentum tensor is T M „ = diag (p, p, p, p) in the comov- 
ing cosmological frame. Therefore, w < — 1 implies viola- 
tion of the NEC. Instability as a consequence of w < — 1 
was studied previously in scalar models 8] and in general 
in Ref. 1]. The analysis applies to dark energy models 
such as k-essence , rolling tachyon condensate ^(j > an d 
phantoms 



II. CLASSICAL FIELD THEORIES 

In this section we will first look at some examples and 
then derive the general result of Ref. which relates 
the violation of the NEC to classical instabilities. Our 
metric convention is such that in a locally inertial frame 
g^v = diag(l, -1,..., -1). 



A. Scalar field model 

Consider the action 

S = J d d x\g\^[Q(X)-V(<P)}, (1) 

where X — <j) a ,^<p a ' fJ ' and <P = (f> a <f) a . Throughout, com- 
mas indicate partial derivatives. By setting the variation 
of the action to zero, we find the equations of motion for 
the fields <p a - 

V^Qx^) = -V^ a . (2) 

In what follows, we assume that we have found a solution 
4> a {x) to these equations of motion, and the quantities 
Qx, etc. are evaluated at this solution. 



1. Stability 

In order to determine whether the solution 4> a (x) to the 
equations of motion is stable against small fluctuations 
S<p a , we consider the second variation of the action, 



8 2 S = 



d d x\g\i (8 2 



8 2 V) 



where 

S 2 Q = (2Q x 8 ab g» v -r^xx 
8 2 V = (2V tt ,S ab + 4V$ <f> a (f> b ) 54> a 



4Q X x<f> a '^ b ' u ) 



(3) 



W4) 
(5) 



We now use a locally inertial frame. For 8 2 C, the 
canonical momentum is 



5jr a = AQ x l 



:a,0 



!xx< 



(6) 



and the Hamiltonian is S 2 H — S 2 K + S 2 V. Here the 
kinetic term is 



8 K = (2Q x 8 ab + AQxx 



a,0 o( Pb,0 

+ (2Qx8 ab 8 v -AQxxftf' 3 ) a 9a<i ocp btj 
where <5</> a ,o are found by solving Eq. |(SJ|, 
1 



8<t>b,j, (7) 



»o,0 



4Qx 

X (5TT b - 



8" - 



2Q X + ^Qxx4>cfi< 



kc.O 



}xx4'b,o4' c ' l 84'c,i) 



(8) 



We are interested in finding a sufficient condition for 
the solutions of this model to be unstable. Moreover, we 
are free to investigate any set of variations we like. In 
particular, let us consider a set of variations such that 
8<Pi,fj, 7^ and 84> a ,n = (a ^ 1). In order to obtain a 
simple condition for instability, we further restrict d — 1 
quantities 8<f>n to satisfy the following equation: 



For this choice of variations, 



8 2 K = 



(9) 



(10) 



We see that the classical solution <j> a {x) is unstable for 
Qx < 0. Specifically, it is unstable to the formation of 
gradients because the system can lower its kinetic energy 
by creating these gradients, while leaving the potential 
energy fixed. See Fig. Q] for a sketch. 



2. NEC 



Now let us consider the NEC for the model JTJ. The 
NEC states that T^n^n" > for any null vector n^. 
The energy-momentum tensor that couples to gravity is 



(11) 




FIG. 1: In a model that violates the NEC, an initially smooth 
field configuration (top) is unstable. The system can lower 
its total energy by developing gradients (bottom), thereby 
decreasing its kinetic energy, while leaving its potential energy 
unchanged. 



Computing the derivative C g ^ = Qx4> a ,u$a,vi we see 
that the NEC requires 



Qx$>^a, M ) 2 >0, 



and therefore, 



Qx >0. 



(12) 



(13) 



Comparing Eqs. (|1(J|) and l|13|) . we conclude that the vio- 
lation of the NEC in this model implies that the classical 
solutions to the equations of motion are unstable to the 
formation of gradients. 



B. Gauge field model 



Consider the action 



S= I d*x\g\* Q(Y)+e\g\-^e^P"F flu F p<7 , 



(14) 



where 9 is a dimensionless parameter, \g\ i^ 9 " is the 
totally antisymmetric tensor in curved spacetime, Y — 
FpuF*", and F M „ = A u . fl - A^ v = A„ t/t - A^ v . A semi- 
colon indicates the covariant derivative appropriate for 
curved space, and the last equality holds because the 
connection has no torsion. Classical electrodynamics cor- 
responds to Q(Y) = -\Y and 9 = 0. 

As before, the equations of motion are found by setting 
the variation of the action to zero and integrating by 
parts. Here the equations of motion for A^ are 



V^Q Y F^ + 9\g\-^e^ a F p , 



0. 



(15) 



In what follows we will assume that we have found a 
solution A^ix) to these equations of motion. 
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1. Stability 



2. NEC 



As in the stability analysis for the scalar field model of 
the previous section, we consider the second variation of 
the action. Since the Lagrangian of (|14|l has no depen- 
dence on the gauge field A^, this is simply 



5 2 S = / d d x\g\i/: Aa ^ A0 JA a;iJl 5Afr v . 



(16) 



Derivatives of C with respect to the covariant derivative 
are equal to those with respect to the regular derivative, 
which in turn are proportional to derivatives with respect 
to the field strength, £ Aa .„ A0 . u = C AailAfiv = \L FaiiFfju . 
Antisymmetry of F afl is used to show the second equality. 
Therefore, to evaluate Eq. (|16[) . we compute the deriva- 
tives of the Lagrangian with respect to the field strength. 
Again using the antisymmetry of F afJi : 



N' 



(17) 



16Q YY F afl F^ 



(18) 
(19) 



where we have defined 

M a/3 = 4Q Y g af3 , 

+89\g\~h a ^. 
The second variation of the Lagrangian now becomes 
5 2 C = (Jlf a V + N 01 ^") 5A a .^5Ap., u . (20) 
For the Lagrangian l|20|) the canonical momentum is 

Sn a = 2M a/3 SA l3 ,o + 2N Mf3h 'SA p . v , (21) 

where we use a locally inertial frame. This allows us to 
write 

S 2 H = (M a(3 + N a °P°) SAa.jaSAp.ja 

+ (M af> 5 l ° -N aif3j )5A a . ti 6Ap. tj , (22) 

and SA a .Q are related to 5ir a and 5A a -i through Eq. (|21Jl . 
To show that this model can be unstable, we choose vari- 
ations such that SAi-fj, ^ and SA a;)J , = (a ^ 1). We 
further restrict d — 1 quantities 8A\-i to satisfy the fol- 
lowing equation: 



N 

Therefore, 



s 2 n 



^SAx-fiSAxfi = N lilj SAi.iSA Uj 



M n SA 1 . SA 1 . + M S l ° SAt-iSAuj 



(23) 



= -4Q> 



(24) 



We see that if Qy > there is the same type of instability 
as in the scalar field case, independent of the value of 8 
in the FF term. Note that in this model the total energy 
is simply the kinetic energy because there is no potential 
energy, i.e., no dependence on the field A^. 



In order to determine whether this model violates the 
NEC, we compute the following derivative 



Cgw — \M a ™ F all Fp v 



F a pF pa , (25) 



where we have used d\g\/dg^ v — — \g\g^ v - From Eq. j2Sfr. 
we see that for this model the NEC requires 



Q Y V a V a > 0, 



(26) 



n^Fa^. Note that the FF term does not 



where W a 

enter into the NEC. We can determine the sign of 'F a 'I ra 
as follows. First notice 



n a <F a = n a n^F a 



0. 



(27) 



by antisymmetry of F a/i . Now use a locally inertial frame, 
and let W a = (%,^) and n a = (|n|,n). Then Eq. 
implies 



n a W a = In^o - |n||#i cos^ = 0, 



(28) 



which means that 'Fa is non-timelike, >F a <P a < 0. There- 
fore, Eq. (|2njl implies that the NEC is violated if Qy > 0. 
As noted below Eq. (|24J) . Qy > is also a sufficient con- 
dition for the model to be unstable. Thus, the violation 
of the NEC in this model implies that the classical so- 
lutions to the equations of motion are unstable to the 
formation of gradients. 



C. General case 

We now give an extended version of the general proof 
of the classical field theory result found in Ref . |lj . Con- 
sider a theory of scalar, 4> a , and gauge, A aa , fields in 
a fixed ci-dimensional spacetime with the metric g plJ . 1 



1 In our analysis we assume a fixed, but possibly non-trivial, space- 
time metric cj MJ/ . We study local instabilities of the matter fields 
propagating on that spacetime background that would manifest 
themselves on microphysical timescales. It is always possible 
to work in an inertial frame at any particular point in space- 
time. One might wonder whether backreaction of the metric 
might dampen the growth of an instability. Since the princi- 
ple of equivalence (which guarantees the existence of an inertial 
frame) applies everywhere, it is easy to see that stability cannot 
be restored unless the null energy condition ceases to be vio- 
lated. That is, although backreaction could in principle decrease 
the magnitude of the negative eigenvalue governing the instabil- 
ity, the eigenvalue cannot go to zero or change sign unless the 
null energy condition is no longer violated. Furthermore, the 
amount by which backreaction might decrease the growth of the 
instability is suppressed by the Planck scale (this is easy to de- 
duce from the Einstein equations), so it is negligible except in 
circumstances in which quantum gravitational effects are impor- 
tant, and known methods cease to apply. 
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We limit ourselves to theories whose equations of mo- 
tion are second order differential equations, so the La- 
grangian density C for the system is assumed to depend 
only on the the scalar fields <j) a , their covariant derivatives 
DpCpa, and the gauge field strengths F afIl/ . The scalars 
may transform in any representation of the gauge group. 
We impose Lorentz invariance on C, but do not require 
overall gauge invariance. That is, we allow for fixed ten- 
sors with gauge indices (but no Lorentz indices) which 
can be contracted with the fields, as is the case, e.g., in 
models with spontaneous symmetry breaking. 

Although we consider minimal coupling between mat- 
ter fields and gravity, non-minimal models of the type 
below are related to minimal models by conformal trans- 
formation Indeed, consider the non- minimally cou- 
pled action 

S = J ' d d x\g\^[\B{ct> a )R + C{<l> a ,D^ a ,F a ^)l (29) 

where R is the scalar curvature. For a single mini- 
mally and non-minimally coupled scalar, B = 1 and 
B = 1 — ^£,<fi 2 , respectively. Under a conformal transfor- 
mation — > = fi 2 g^, the action l|29l) transforms 
as follows: 

S^S = J d d x\g\?Q d {£ + \Bf2- 2 [R 

-2(d- l)DlnJ?- (d-2){d- l)(Vln/2) 2 ]},(30) 

where C is within the same class as C, and all indices 
are raised and lowered with the untransformed metric. 
Notice that if we choose fl = f^A 2 -"*) and integrate Eq. 
(I30J1 by parts, the conformally transformed action be- 
comes the action for minimally coupled scalars with a La- 
grangian density of no larger generality than C. (Deriva- 
tives of B give terms involving gradients of the scalar 
fields; after adding these with additional multiplication 
by a -dependent terms, the resulting Lagrangian is still 
within a class of C considered here.) The conformally 
transformed action is 

S = J d d x\g\i[^R + B^£-^^(VlnB) 2 }. (31) 

Thus, there is no need to consider non-minimally coupled 
models of Eq. (|2*9")l separately. 
For the action 

S = J d d x\g\?£((b a ,D^ a ,F a ^) (32) 

to be stationary, its first variation has to vanish, 5S = 0. 
This leads to the equations of motion for the fields 4> a 
and A aa ; in the classical analysis we assume that we 
have found solutions to these equations, about which we 
expand. See Fig. |2Jfor a flow chart of the method that we 
use to show a connection between violation of the NEC 
and instability. 





FIG. 2: A flow chart of the general version of the proof that for 
models described by the action (1321 . only solutions satisfying 
the NEC can be stable. 



1. NEC 



We use the notation ipA^ = (D^<f) a , F aa ^), where the 
abstract index A runs over both Lorentz and color in- 
dices, as well as type of field. For example, in Sec. Ill Al 
the index A is a color index, ipA^ — <t>a,n, while in Sec. 
Iflfil it is a Lorentz index, ipAu = Fan- I n general, for 
a model with A^ s scalar fields and N e gauge fields in d 
dimensions, A runs from 1 to N s + |(d- 1)N S . 

Assuming that ipA^ and g^ v are independent variables, 
we now prove that there is a relation between the deriva- 
tives of C with respect to them. Since we consider only 
Lorentz invariant models, the Lagrangian can be written 
in terms of a number of Lorentz invariant quantities: 



C = C(X U X 2 ,...). (33) 



Objects represented in Fig. [3] are the simplest examples 
of X n . Each dot represents a Lorentz index, and a line 
connecting them denotes contraction using the metric. 
Rectangles (with two indices) are field strengths, small 
circles covariant derivatives of scalar fields, and large cir- 
cles epsilon tensors. All Lorentz indices are ultimately 
contracted, and we suppress color indices for simplicity. 
As we will show, the relation between the derivatives of 




FIG. 3: Representation of simple Lorentz invariants that the 
Lagrangian may depend upon. Each dot represents a Lorentz 
index, and a line connecting them denotes contraction using 
the metric. Rectangles (with two indices) are field strengths, 
small circles covariant derivatives of scalar fields, and large 
circles epsilon tensors. All Lorentz indices are ultimately con- 
tracted, and we suppress color indices for simplicity. 



C is 



n 

n 

\M AB i> A ^ Bv + y^K, (34) 



(M [n) AB ^ B ^ + ^- v *L (n) J 

n 



(35) 



where quantities Mr n \, Kr n ) and Li n ) originate from the 
dependence of £ on X n . The first term in Eq. (|3"4*|) is ob- 
tained by noting that for each and every g^ v in C there 
are two -i/^s attached to it. The second term is present 
if C has a term involving the curved space totally an- 
tisymmetric tensor \g\~^e Vl '" Vd , again recalling that the 
derivative of the determinant of the metric is propor- 
tional to the metric. Similarly, differentiation with re- 
spect to ipAn yields the M and L terms in Eq. 135(1 . 

To prove this relation for the most general case, we 
consider the most general Lagrangian of the type consid- 
ered in this paper, represented by Fig. 21 The notation is 
the same as in Fig. with the addition that Z denotes 
the elements of the Lagrangian which are shaded gray. 
We focus our attention on the term represented by the 
elements which are shaded black. 

Indeed, consider the portion of Fig. 0] whose elements 




1 1 1 1 1 1 1 



FIG. 4: Representation of the most general Lagrangian of 
the type considered in this paper. The notation is the same 
as in Fig. [3] We focus on the elements shaded black. In 
graphical terms, M for this term in the Lagrangian can be 
obtained from the figure by simply removing these elements. 
Z represents the remainder of the Lagrangian. 



are shaded black; it equals 



C=\g\-h a -F atl g^F v(3 g^Z 1 



(36) 



(Again, we suppress color indices for simplicity, as they 
do not affect the proof.) By differentiation, we obtain 
Eqs. (gll and with 



M aP = -\g\-? (e— sr> + f?~sr<) Z-y-> (37) 
K = \g\-ic a -F ap gP-F a0 g^Z T .., (38) 
L a ... = -\g\-ig ap F P0 g^Z T ... (39) 

Note the derivative Cf u « generates a contribution to M 
which is matched by a corresponding contribution from 
C g 0-y. Other contractions of fields with g^ v (i.e., as indi- 
cated in the figure) can be analyzed similarly. In graph- 
ical terms, M can be obtained from the figure by simply 
removing the elements shaded black. 
From Eq. the NEC requires 



V A M AB 1> B > 0, 



(40) 



where <F A = ipA^. Thus, to satisfy the NEC, M AB has 
to be positive semidefinite on the subspace spanned by 
This property is crucial for stability of solutions, to 
which we now turn. 



2. Stability 

To study the stability of the solution tpA (x) , we con- 
sider the second variation of the action, 



S S = J d x\g\* {C^ A ^ B 8ip A 8ip B 
+ 2C^ a . 4 , b x 8iJ)aHb:\ + C^a^b^Ha-uHb-v) ■ (41) 
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Here quantities £ 



dC/dipA, etc. are evaluated at 



iPa(x)- Also notice that -0A ;M = (D M tj) a , A ao 
covariant derivatives of ipA, are different from i^a^ 

Differentiating Eq. I|35() we find 

M AB g^ + N A ^ By . 



the 



(42) 



The separation of the second derivative into the first and 
second terms in Eq. i|42|) is natural: g^ v appears only 
in the first term, and iV represents all remaining terms. 
N contains terms obtained by differentiating M and L 
with respect to ipBv, plus additional terms if -0 is a field 
strength, as in Eq. H19JI • The v index obtained from these 
ij)Bv derivatives is attached to a field and not the metric 
g^ v . (L does not contain an epsilon tensor since Z does 
not.) 

Next notice that even though ipAp, and ipAm differ, the 
derivatives of C with respect to them are related due to 
the form of the action l|32[l . Thus, the second variation 
of the Lagrangian in Eq. (|41l) becomes 

S 2 C = (M AB g»» + N a ^ Bv )8^ A ;Mb- <v 

+ 2£ Mb .J*1> a 5^b;X ~ S 2 V, (43) 

where 5 2 V = -C^, a ^ b 8^ a 5^b- 

Let us use a locally inertial frame. For the La- 
grangian (|43|l . the canonical momentum is 

6n A = 2M AB 5^ B , Q + 2N A0Bv 5^ B -u 

+2L MA;0 5ip B , (44) 

which leads to the following effective Hamiltonian for 
fluctuations about the classical solution: 

S 2 H = S 2 K + S 2 V - 2L MB ..5il) A WB;i> (45) 

where the kinetic term is 

8 2 K = (M AB + N A0B0 ) 5^A; S^B- fi 

+ (M AB 6 l > - N AiB i) 6Tp A ;iSlpB;j. (46) 

Here <50A;O are functions of Sitb, S^b, and SipB-,i as found 
from Eq. JUJl. 

We now prove that nonnegativeness of the kinetic term 
S 2 K implies positive semidefiniteness of the matrix M AB . 
Indeed, suppose that M AB is not positive semidefinite. 
In such a case, the matrix M AB has at least one negative 
eigenvalue, which means that there is a basis in which 
the matrix is diagonal with at least one negative entry, 
M AB = diag (mi, . . . , m„) (mi < 0). Let us choose field 
variations which are nonzero only in the direction of the 
negative eigenvalue: <50i ;A1 ^ and StpA-^ = (A > 1). 
We further restrict d — 1 quantities 8^>\-i to satisfy the 
following equation: 

7V lolo <50i ;O <5Vi;o = N ul '5^,Mi;o- (47) 

These conditions make the kinetic term in Eq. i|45|) neg- 
ative, 



S 2 K = m : 



Wl;0) 2 



thus proving that in order for S 2 K to be nonnegative, the 
matrix M AB has to be positive semidefinite. 

Note that if the background field leaves an unbroken 
rotational symmetry — i.e., there are d — 2 or fewer inde- 
pendent vacuum expectation values if)^ — one can find a 
solution to Eq. (|T7|) in any frame, by taking <50i ; o = 
and the spacelike fluctuations to have support only in the 
direction(s) orthogonal to all of the ip^- hi the remaining 
cases (which, in particular, have no unbroken rotational 
symmetries) there may be Lorentz frames in which no so- 
lution to Eq. ljT7ll exists. The examples given by Ref. @ 
are in this category. In Refs. and 0, it is emphasized 
that superluminal modes play a crucial role in permit- 
ting stability and NEC violation simultaneously. How- 
ever, models with superluminal modes are not causal. In 
the appendix, we show that the assumption of causality 
is strong enough to restore the direct connection between 
NEC violation and instability. 

If the kinetic energy is negative, then the system de- 
scribed by the Hamiltonian (|45|l is (locally) unstable. 
The linear term in Eq. I|45|) can never act to stabilize 
the system. Therefore, in order to determine whether 
the system is unstable, we focus only on the quadratic 
terms. If the potential energy is positive then small per- 
turbations will cause the classical solutions to grow expo- 
nentially away from the original stationary point. How- 
ever, it is possible to have classical stability if one chooses 
the potential energy to be negative; in this case we have 
an upside-down potential with negative kinetic term, or 
a "phantom" . Such models necessarily exhibit quantum 
instabilities ^l| . Using the result established in the pre- 
vious paragraph, we thus conclude that solutions to the 
theory given by the action (|32|1 are stable only if the ma- 
trix M is positive semidefinite. 

Combining the relations between nonnegativeness of 
8 2 K and positive semidefiniteness of M AB on one hand, 
and the NEC and positive semidefiniteness of M AB on 
the the other hand, we conclude that for the theory given 
by the action i|32|) , only solutions satisfying the NEC can 
be stable. 



III. QUANTUM FIELD THEORIES 

We can deduce similar results for quantum systems. 
Suppose there exists a quantum state \a) and a null vec- 
tor ?i M such that the NEC is violated in a quantum aver- 
aged sense: 



(alT^laJn'W = (a\M AB <P A V B \a) < 0. (49) 

Define a basis \cf>) in which the operator M. = M AB] I> a^b 
is diagonal: 



< 0, 



(48) 



M\4>)=m(<p)\<p). 



(50) 
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By orthonormality of |</>), we see that violation of the 
quantum averaged NEC implies 

(a\M\a) = ^(aMWMWMW 
W 

= ]T I(«|0)| 2 ™(0) <0. (51) 



This means that there exist eigenstates \<j>), whose over- 
lap with | a) is non-zero and on which the operator M. has 
negative eigenvalues. This requires that M AB and hence 
8 2 K is not positive semidefinite; by continuity, this must 
also be the case in a ball B in the Hilbert space of \<j>). 

As a further consequence, we can conclude that a state 
|a) in which the quantum averaged NEC is violated can- 
not be the ground state [l^- Suppose that |a) is an 
energy eigenstate: H\a) = E a \a). An elementary result 
from quantum mechanics is that |a) can be the ground 
state only if 



E a = (a\H\a) < (a'\H\a') 



(52) 



for all normalized states \a!) which need not be energy 
eigenstates. However, it is possible to reduce the expecta- 
tion value of H by perturbing \a). Specifically, we adjust 
| a) only in the ball B, where we know from Eqs. (14611 1481) 
that there are perturbations which reduce the expecta- 
tion of the kinetic energy without changing the expec- 
tation of the potential. For a pictorial representation of 
this, see Fig. ^ where now the two different graphs are 
sketches of two eigenstates of the operator M.. Let the 
top one be \4>\) and the bottom \<j>2)- Expand the state 
\a) in the basis of these eigenstates: 



!«> = £< 



(53) 



Since by assumption the state \a) violates the quantum 
averaged NEC, the magnitude of coefficients of smooth 
eigenstates, e.g. \4>i), will decrease, while those of eigen- 
states with large gradients, e.g. \<j>2), will increase, until 
the quantum averaged NEC is no longer violated. Thus, 
a state \a) which violates the quantum averaged NEC 
cannot be the ground state because perturbations within 
the ball B can lower the expectation value of the energy. 

Note that the discussion above is in terms of unrenor- 
malized (bare) quantities. The renormalized expectation 



(a\M mn \a) = (a\M\a) - (0\M\0), 



(54) 



where |0) is the flat space QFT ground state, could be 
negative (e.g., as in the Casimir effect |2j), but this is 
possible only if |a) is not the ground state. 

In known cases of NEC violation, such as the Casimir 
vacuum or black hole spacetime, it is only the renormal- 
ized energy- momentum tensor which violates the NEC. 
As a simple example, consider a real scalar field cf>. The 
energy-momentum tensor is simply T^ v = d^<j)d u (j) plus 
terms proportional to g M „ which do not play a role in the 



NEC. Then M. = (n^d^4>) 2 is a Hermitian operator with 
positive eigenvalues. Therefore, its expectation value in 
any state is positive: 



(a\M\a) > 0, 



(55) 



for any a), including the Hartle-Hawking, Casimir or 
flat-space vacuum. We can verify this by direct cal- 
culation in d = 4 dimensions, computing the energy- 
momentum tensor using point-splitting regularization: 



(0\<P(x)<p(x')\0) 



1 



4tt 2 | 



which yields 



(0\T ltu (x ) a/)\0)nf l n v = 



2[n^(x - x')^ 

TT 2 \x — x'\ 6 



(56) 



(57) 



So (0|jM|0) is manifestly positive. Note that this result 
holds even if the expectation is taken in an arbitrary 
state rather than |0), since the bare expectation is always 
dominated by the UV contribution. Now, had we taken 
a negative kinetic energy term for the scalar, the overall 
sign of M. changes, allowing violation of the NEC. But, 
this model is clearly unstable, in accordance with our 
results. 

Our results are summarized as follows, where \a) is any 
quantum state, and only models of the type described by 
action are considered: (1) (alA^a) < 0, i.e., the 
bare NEC is violated, implies that \a) is locally unsta- 
ble, (2) (alA'Jrenla} < 0, i.e., the renormalized NEC is 
violated, implies that \a) cannot be the global ground 
state. 



IV. MODELS WITH FERMIONS 

In this section we present a generalized version of 
the treatment of fermions in Q. Both Dirac and Weyl 
fermions with standard couplings to the bosonic fields 
of Sec. IIIII may be included in the model without any 
effect on the results of that section. In the purely classi- 
cal analysis, fermions play no role, while in the quantum 
case we can integrate them out to obtain their contribu- 
tion to the energy-momentum tensor. We will see that 
this fermionic contribution always satisfies the NEC. 

In addition to the bosonic part in the Lagrangian, we 
include the following fermionic part: 



£« =$[i0-m(^ a )]V. 



(58) 



A Dirac fermion tp couples to scalar and vector bosons 
via a covariant derivative = + and an arbitrary 
real scalar function m(<f> a ). We use Euclidean quantities 



obtained via substitutions x 

-r 



o 



-ix°, A -> iA , 7° 



DtpDifj exp 
det (i Tp — m). 



- / d d x\g\^C {i) 



(59) 



8 



The Euclidean operator ft is self-adjoint. 

The imaginary part of the action vanishes since it is 
an odd function of ft and the trace of odd number of 
gamma matrices vanishes: 

Imf = -±tr In [(i ft - m)(-i ft - m)^ 1 ] 

= itrarctan (ftm^ 1 ) = 0. (60) 

Using eigenvectors ip k and eigenvalues A^ of the operator 
i ft — m, the real part of the action becomes 



ReF = -±tr \n[(-i ft - m)(i ft - m)\ 



(61) 



Since the operator i ft — m is not self-adjoint, the eigen- 
values are complex and the eigenfunctions are not orthog- 
onal. However, they can be normalized, J d d x\g\^i/j k ijj k = 
1. Using this normalization we have 



|A fe | 2 = 



x\g 



hli-ift 



m)(i ft - m)i>k. (62) 



For the purposes of calculating the contribution of Re F 
to the energy momentum tensor, it is enough to consider 
only those terms in Eq. l|62|l which contain g^ v : 



I A* 



d d x\g\ig^(D^ k )\D^, 



(63) 



We have used here anti-hermiticity of £)„• Thus, using 
the definition of the energy momentum tensor 



t!S = 2 



- fiiy 
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. ST 

JgT^ 1 



(64) 



the fermionic contribution to the NEC expression is non- 
negative: 



k 



\\ k \- 2 (n»D^ k )\n v D^ k )>Q. (65) 



The Weyl operator does not define an eigenvalue prob- 
lem, as it maps a spinor into the opposite chirality; see 
e.g. However, the Weyl determinant can be defined 
as the square root of the corresponding Dirac determi- 
nant obtained by considering the Weyl fermion and a 
partner in the complex conjugate representation. The 
real part of the contribution to the energy-momentum 
tensor is exactly one half of the Dirac contribution con- 
sidered previously. The presence of either type of fermion 
does not lead to violation of the NEC, and so the results 
of Sec. IIIII hold in models with arbitrary fermionic and 
bosonic particle content. 



V. FLUIDS 
A. Perfect Fluids 

A macroscopic system may be approximately described 
as a perfect fluid if the mean free path of its components 



is small compared to the length scale of interest. For the 
dark energy, this length scale is of cosmological size. A 
perfect fluid is described by the energy-momentum tensor 



(66) 



where p and p are the energy density and pressure of 
the fluid in its rest frame, and is its velocity. Let 
= Ju^ be the conserved current vector (j^.^ = 0), 

and J = (ijuj 71 ) 5 the particle density. 

The energy-momentum can be written 0, as 



(67) 



where, comparing with Eq. I)66|l . we have p = f(J) and 
p = Jf — f. The function /(J) implicitly determines the 
equation of state. Also, prime indicates derivative with 
respect to J at fixed temperature and volume. 
The NEC for the tensor (JHZJl becomes 



T^nW = (f/J)(j^) 2 > 0. 



(68) 



Thus, perfect fluids with negative /'(J) violate the NEC. 
Below, we demonstrate that /'(J) < implies an insta- 
bility. 

Recall that the speed of sound in a fluid is given by 
c s = (dp/ dp) i — (J/"//0 5 ) an d that complex c s (or 
neg ative compressibility) implies a mechanical instability 
jl5j . This instability is avoided only if /" < 0. However, 
if /' and /" are negative, then the fluid is unstable with 
respect to clumping. To see this, we first recall that 
the Helmholtz free energy F is a Legendre transform of 
the internal energy E, i.e., F = E — TS, where S is the 
entropy. In the thermodynamic limit, entropy scales with 
particle number N = JV, and so 



F = Vf - JVTs, 



(69) 



where from above, f(J) = p is the energy density, and 
s = S/N is the specific entropy, which is a function of 
temperature but not particle number. Therefore, 



d 2 F 



dJ 2 



(70) 



T.V 



Now consider two adjacent regions of the fluid with 
identical volumes. Suppose we transfer a small amount 
of matter 5 J from one volume to another. Expand the 
change in free energy to second order in the number den- 
sity: 



6 Fx + SF 2 



dF 
dTx 
dF 
~dFi 



SJx 



SJ 2 



ld 2 F 
2~dJl 
ld 2 F 

2~ojJ 



(SJx) 2 



(SJ2) 2 



(71) 



The linear terms cancel due to particle number conser- 
vation, SJx — —5J 2 — 5J. Therefore, using Eq. lj?0"]l. the 
resulting change in total free energy is given by 



5F = ±Vf"(J)(5J) 2 < 0. 



(72) 
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We see that the system can decrease its free energy by 
clumping into over- and under-dense regions. This itself 
is an instability, which results in a runaway to infinitely 
negative free energy unless the assumption of negative /' 
or negative /" ceases to hold. NEC violation cither leads 
complex speed of sound or clumping instability. 

B. Beyond Fluids 

We can generalize our analysis beyond fluids, relaxing 
the assumption of the existence of a conserved current 
j M , and requiring only energy conservation and statistical 
equilibrium. 

Consider an isolated system at fixed energy E. If the 
system is ergodic it can (after a sufficiently long time) be 
described by a microcanonical ensemble, which uniformly 
samples phase space. If it has many degrees of freedom, 
the microcanonical ensemble is well approximated by a 
canonical ensemble with a Boltzmann weighting, in which 
the energy is allowed to fluctuate. In this case the ther- 
modynamical relation E = TS — pV is satisfied, with 
T = (dE / dS)\v ■ Dividing by the volume, we obtain 
p + p = Ta, where a = S/V is the positive-definite en- 
tropy density. Thus, for a system which is homogeneous 
and isotropic, i.e., one with energy- momentum tensor as 
in Eq. I|66(l . violation of the NEC requires negative tem- 
perature, or entropy (number of states) which decreases 
with energy. For a field theory, this suggests either a 
negative kinetic energy term or an "upside-down" poten- 
tial; if the kinetic and potential energies have opposite 
sign (the usual case), the number of states available to 
the system increases as energy is added, and the tem- 
perature is positive. Moreover, in order for the partition 
function (trace over Boltzmann weights) to converge, the 
kinetic energy term must be negative if the temperature 
is. This is because at any fixed potential energy there is 
an infinite volume of phase space with arbitrarily large 
gradients, i.e., kinetic energy. Therefore, if the kinetic 
energy term and the temperature are not of the same 
sign, the partition function does not converge. 

An example of a condensed matter system with a neg- 
ative temperature is a ferromagnet in an external mag- 
netic field [l6(. Beyond a certain threshold energy, the 
number of states available to the system decreases as en- 
ergy is added, i.e., the temperature is negative. Due to 
the presence of an external magnetic field, however, the 
system is anisotropic; its energy-momentum tensor is not 
of the form of Eq. 1|66[) . and there is no connection be- 
tween negative temperature and violation of the NEC. It 
is only for homogeneous, isotropic systems, e.g., systems 
of interest in cosmology, that this relation exists. 

VI. CONCLUSION 

We have exhibited a direct connection between viola- 
tion of the NEC and instability in a wide variety of mod- 



els. One important application is to the dark energy, 
whose equation of state w should be no more negative 
than —1 in order not to violate the NEC. 

In classical field theory — including all Lorentz- 
invariant models involving both minimally and non- 
minimally coupled scalar and gauge fields with second- 
order equations of motion — violation of the NEC gen- 
erally implies the existence of an instability toward the 
formation of gradients. In the case of configurations that 
possess at least some rotational symmetry (e.g., those of 
interest in cosmology), violation of the NEC always im- 
plies instability (Sec. IIC). In cases with less symmetry, 
it is possible for NEC-violating models to be stable if 
they exhibit superluminal excitations |(|; however, such 
models are not causal Q. In the appendix we show that 
in any causal scalar model, NEC violation implies insta- 
bility. 

These results can be generalized to quantum theories, 
including fermions. Analogous results for perfect fluids 
relate violation of the NEC to either a complex speed of 
sound or a clumping instability. 
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APPENDIX: CAUSALITY AND STABILITY 

In this appendix we impose an additional causality 
constraint, which excludes superluminal excitations of 
the type considered in Refs. @, Q. We show that a 
causal, NEC-violating model cannot be stable. For sim- 
plicity, we limit our attention to the case of scalar fields. 

1. Fluctuations 

Consider a theory of n scalar fields <f> a in d-dimensional 
spacetime. Using the first derivative of the field, <p afl = 
V M </) a , we can construct a Lorentz-invariant matrix 
X a b = '4> a ^4>bu- The most general Lorentz invariant 
Lagrangian resulting in the equation of motion of at most 
second order is C = F(X a i > , cf) c ), where the function F is 
a Lorentz scalar, but in general it is not a scalar in the 
flavor space. The second fluctuation of the Lagrangian is 

6 2 C = F a ^S4>aMbu + ■ ■ ■ , (A.l) 

where 

F a, bv = Fx ^u + 2F XacXbi ^ d . (A.2) 
The equation of motion for the fluctuations is 

F a ^5cf> b ^ + • • • = 0. (A.3) 
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In Eq. (|A.1(1 . ellipses denote terms at most linear in the 
derivatives of fluctuations; such terms are not important 
for the present analysis. Specifically, they do not enter 
into the equation for the characteristic surfaces or the dis- 
persion relation for short wavelength fluctuations. See, 
e.g., Ref. [13 for more on characteristics in partial differ- 
ential equations. In the case of the dispersion relation, 
the short wavelength approximation dictates that only 
terms with the largest number of derivatives matter. 

2. Hyperbolicity and causality 

When the equation of motion is hyperbolic, its do- 
main of influence is bounded by the outside character- 
istics originating from the region where nontrivial initial 
data is prescribed. The resulting characteristic surface is 
given by £(x) = 0, where £ M = V M £ satisfies the charac- 
teristic equation 

11^^11=0, (A.4) 

and ||--- || denotes a determinant in flavor space. For 
a hyperbolic model, characteristics and therefore £ M are 
real. For such a model, the characteristic surface is inside 
the light cone £q — < (i.e., the model is causal), 

if and only if the matrix F a ^ bu n^n v is sign-definite for all 
null n^. By writing F a » bv = J2 C F c Vv >^ wh ere Fg u and 
u" are eigenvalues and eigenvectors, respectively, this re- 
sult is obtained by noticing that the characteristic surface 
consists of the null surfaces for the metrics Fg v . Each of 
these null surfaces has to be inside the light cone, hence 
the above result. From Eq. (|A.2(1 . this means that the 
matrix C ab — 2Fx ac x bd {4>c n v)( < l ) d n v) ^ s sign-definite. 

In terms of A ab = F Xab , B ab = 2F Xac x bd (€W d ^), 
and w = —^n^, the characteristic equation is 

\\B-luA\\=0. (A.5) 

For nonsingular A, the parameter u> is equal to one of 
the eigenvalues of the matrix A~ 1 B. Causality requires 
oj > and so 2 A~ X B > 0. 
Consider a Lorentz boost 

» - « kn„ (A.6) 

for some particular null vector n^. (It is easy to check 
that £ M can always be made nearly null by appropriate 
choice of frame.) In this new frame, B ab is arbitrar- 
ily close to the matrix C ab evaluated on this particular 
rip. Then, because C is sign-definite for all n M , the sign- 
definitcness of B follows. This means that if a model is 
causal, then either A > 0, B > or A < 0, B < 0. This 
result can be extended to singular A as well. 



2 The inequality Q > (Q < 0) for the matrix Q means that Q is 
positive (negative) definite. 
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FIG. 5: The NEC, causality (here C), and quantum mechani- 
cal stability (the absence of ghosts, here (ft) are all determined 
by the sign-definiteness of matrices A and B. The NEC holds 
for A > 0; the model is causal for A" 1 B > and ghost-free 
for B > 0. A model is ghost-free and causal only if it obeys 
the NEC. 



3. Stability and NEC 

To study local stability, it is enough to consider the 
equation of motion in a small neighborhood of a given 
point x. In such a neighborhood, we can ignore spacetime 
dependence of F a>lbv and solve the equation in terms of 
Fourier modes with momenta k^. Also, in the short wave- 
length approximation, only terms quadratic in deriva- 
tives matter. The corresponding dispersion relation is 

llF^fc^ll =0. (A.7) 

Classical stability requires that /c M is real, which means 
that the matrix F afibl ' satisfies the same condition which 
is required for the characteristics to be real. In other 
words, models are classically stable if and only if they 
are hyperbolic. To avoid ghosts, i.e, quantum mechani- 
cal instabilities , there is an additional condition that 
the matrix F a0b0 be positive definite. Except in the triv- 
ial case where Fx ab x cd vanishes, we can always choose 
such a Lorentz boost that F a0b0 is arbitrarily close to 
(jab £ or some n ^ Since B and C are of the same sign 
definiteness, the matrix F a0bQ can be positive definite, 
and therefore the model is ghost-free, only if B > 0. 
Then, from the previous section, the combined require- 
ment that the model be both causal and ghost-free can 
be satisfied only if A > and B > 0. 
From the energy-momentum tensor 

T MV = -Fg^ + 2F Xab <t>a»<t>b»> (A.8) 

it follows that the NEC is equivalent to A > 0. Thus, for 
a hyperbolic, causal (A^B > 0), classically and quan- 
tum mechanically stable model (B > 0), the NEC is sat- 
isfied (A > 0). Vice versa, a causal, hyperbolic (and 
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therefore classically stable), NEC- violating model can- 
not be quantum mechanically stable. See Fig. [S] for a 
diagram summarizing the relation between quantum me- 
chanical stability, causality, and the NEC. 



the absence of ghost instability imply non-causal behav- 
ior for this model. Moreover, for Fx < and Fxx > it 
is impossible to find a Lorentz frame in which the Hamil- 
tonian 



4. Example 

For the case n = 1, combining Eqs. (| A. 2|> and (|A.4(1 . 

the characteristic equation is 



s 2 n 



F x + 2Fxx(rr Wo) 



F X S» - 2F X x4> t fi> 



(A.10) 



FxH^ + ZFxxir^) 2 = 0, 



(A.9) 



and the model is causal if F^Fxx > 0. Violation of the 
NEC and absence of ghost instability imply Fx < and 
Fxx > 0, respectively. We see that NEC violation and 



is positive definite since the second term is always neg- 
ative. Clearly, this can happen only for a non-causal 
model. When superluminal modes are present (i.e., a 
model is not causal), the sign of the Hamiltonian can be 
changed by performing a Lorentz boost. 
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